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An exact solution of two singularities in the teleparallel equivalent to general relativity theory 
has been obtained. A holographic visualization of the binary black holes (BBHs) space-time, due to 
the non vanishing torsion scalar field, has been given. The acceleration tensor of BBHs space-time 
has been calculated. The results identify the repulsive gravity zones of the BBHs field. The total 
conserved quantities of the BBHs has been evaluated. Possible gravitational radiation emission by 
the system has been calculated without assuming a weak field initial data. 
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I. INTRODUCTION 

The discovery of the first pulsars in a binary system 
PSR B1913-I-16 [1] represents a natural laboratory for 
testing theories of gravity. As a matter of fact the param¬ 
eterized post Newtonian technique recommends the GR 
theory to pass the binary pulsar laboratory test. How¬ 
ever, the general relativistic exact solution of such sys¬ 
tems is not recognized yet. The first trial was by [2] as 
an axially symmetric vacuum solution of two monopoles, 
the solution has been reobtained in [3]. A more exten¬ 
sive studies on the two body problem in GR has been 
studied [4-9]. But these solutions, in fact, are in a series 
trouble with GR theory. The so-called Gurzon solution 
has been obtained as a vacuum solution of two singular¬ 
ities, while the solution violates the elementary flatness 
condition between the two monopoles. The only physical 
interpretation is that the region between the singulari¬ 
ties cannot represent a vacuum solution of Einstein field 
equations by introducing an inert strut (i.e conical sin¬ 
gularity) to keep the static configuration of the system 
balanced [10, 11], while other trials to remove the strut 
was by considering the spin-spin interaction [12], later it 
has been shown that double Kerr space-time cannot af¬ 
ford a balanced system [13]. Also, the solution shows a 
peculiar singularity behaviour, whereas it exhibits a di¬ 
rectional singularity allows the particles approaching the 
singularity along the axis of symmetry to be geodesically 
complete accessing to some new region of the space-time 
[14]. This leads some to conclude that the singularities in 
Gurzon solution are actually rings rather point-like par¬ 
ticles [15-17]. Interestingly, for the superposition of a 
ring and another body, a membrane-like singularity ap¬ 
pears besides the strut one [18, 19]. Also, it has been 
shown that the interior of cylindrically symmetric space- 
times topology is physically acceptable, if there exists 
a negative energy cosmic string density [20]. Moreover, 
this solution has been modified to describe the pulsation 
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time delay (Shapiro delay) of signals from binary pulsar 
systems by using a co-rotating coordinates [21]. 

In our study, we investigate possible solution of two 
monopoles in the TEGR framework. We expect new in¬ 
sights and interpretations when use different tools, i.e 
teleparallel torsion gravity, of such systems. One of the 
important features of the TEGR theory, that is the vac¬ 
uum solutions do not imply the vanishing of the telepar¬ 
allel torsion scalar, unlike the GR case in which the van¬ 
ishing of the Ricci scalar in vacuum solutions is necessar¬ 
ily. So we expect the torsion scalar to play an important 
role to visualize the space-time, the more important role 
is to provide an alternative source of the compressional 
medium between the two singularities. We organize the 
article as follows: In Section II, we briefly review the La- 
grangian formalism of TEGR theory. In Section III, we 
derive an exact solution of the TEGR field equations for 
a BBHs space-time. In Section IV, we use the non van¬ 
ishing torsion scalar to visualize the vacuum solution of 
the BBHs gravitational field. In Section V, we evaluate 
the acceleration tensor to determine the repulsive gravity 
zones of the BBHs field. In Section VI, we calculate the 
conservative quantities of the BBHs space-time and an 
estimate value of its gravitational radiation emission. In 
Section VH, we summarize and conclude the work. 


II. TELEPARALLEL EQUIVALENT OF 
GENERAL RELATIVITY 

The TEGR is an alternative description of gravi¬ 
tation of Einstein’s GR. The theory has been used 
to reexamine many GR solutions [22-24]. It is con¬ 
structed from tetrad helds^ (e'^i) instead of a met¬ 
ric tensor fields The metric tensor can be con¬ 
structed from the tetrad fields: gij = where 

O/ii, = diag(I, — I, — I, — I) is the Minkowski metric for 


^ The Greek indices describe the components of tangent space to 
the manifold (space-time), while the Latin ones describe the com¬ 
ponents of the space-time 
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the tangent space, so the Levi-Civita (symmetric) con¬ 
nection can be constructed. However, it is possi¬ 
ble to construct Weitzenbbck nonsymmetric connection 
[25]. The Weitzeinbock 4-space is de¬ 
scribed as a pair (M, e^), where M is an 4-dimensional 
smooth manifold and (/r = 0, • • • , 3) are 4-linearly 
independent vector fields defined globally on M. The 
Weitzenbbck space is characterized by the vanishing of 
the tetrad’s covariant derivative, i.e. V.e^i = 0, where 
the covariant derivative is with respect to (w.r.t.) 
the Weitzenbbck connection. So this property identifies 
auto parallelism or absolute parallelism condition. As a 
matter of fact, the operator is not covariant under 
S'0(l,3) group, i.e local Lorentz transformations (LLT). 
The lack of the local Lorentz symmetry allows all LLT 
invariant geometrical quantities to rotate freely in every 
point of the space [26]. In this sense, the symmetric met¬ 
ric (10 degrees of freedom) cannot predict exactly one 
set of tetrad fields; then the extra six degrees of freedom 
of the 16 tetrad fields need to be fixed in order to iden¬ 
tify exactly one physical frame. It can be shown that 
the metricity condition is fulfilled as a consequence of 
the absolute parallelism condition. The Weitzenbbck is 
curvature free while it has a non-vanishing torsion given 

by 


rpl . -p2 -pZ 

-L jk ’— i jk -L kj • 

(1) 

and contortion is given by 



(2) 

In the teleparallel space one may define three 
Weitzenbbck invariants: Ii = I 2 = 


and /a = T'T^, where T* = Tjb. We next define the in¬ 
variant T = All + BI 2 -f C/a, where A, B and C are ar¬ 
bitrary constants [26]. For the values: A=l/4, 5 = 1/2 
and C = —\ the invariant T is just the Ricci scalar , 
up to a total derivative term as we will show below; then 
a teleparallel version of gravity equivalent to GR can be 
achieved. The teleparallel torsion scalar is given in the 
compact form 

(3) 

where the superpotential tensor is defined as 

:= i + 6iT^\ - , (4) 

which is skew symmetric in the last two indices. We next 
highlight some useful relations between Riemannian and 
teleparallel geometries. The double contraction of the 
first Bianchi identity of the teleparallel geometry gives 

=-TT) _2V,rA^-, (5) 

where the covariant derivative V is w.r.t. the Levi-Civita 
connection T. The second term in the right hand side 


is a total derivative. So, it has no contribution in the 
variation when use the right hand side instead of the 
Ricci scalar in the Einstein-Hilbert action. Consequently 
the Ricci and teleparallel torsion scalars are equivalent up 
to a total derivative term. In spite of this quantitative 
equivalence they are qualitatively different. For example, 
the Ricci scalar is invariant under LLT while the total 
derivative term is not, so the torsion scalar. Accordingly, 
the TEGR Lagrangian is not invariant under LLT [27- 
29]. 

The TEGR gives the action a gauge gravitational field 
Lagrangian as [30, 31] 

/Vf2 r 

S=^ J \e\[T + £Mi<i>A)] ( 6 ) 

where Cm is the Lagrangian of the matter fields ‘I’a and 
Mpi is the reduced Planck mass, which is related to the 
gravitational constant G by Mpi = yJhcj'inG. Assume 
the units in which G = c = h = 1. In the above equation, 
|c| = '\/~9 = det(e^i). The variation of (6) w.r.t. the 
tetrad fields e^i give rise to the following field equations 
[30] 

9,(eV')=4^ee^^(t/ + 0/), (7) 

where = e^J'SCC the (pseudo) tensor tA is 

= ^{^T^rmSl^^-5lT), (8) 

and the matter energy-momentum tensor is 

As the tensor 5 ^^ jg skew-symmetric, i.e 
this implies that didj{eS^'‘^) = 0 [26]. Therefore, 

dj [eefj,\ti^ + ©i-’)] = 0. 

The pseudo tensor tj has no equivalent in GR theory. To 
reveal its nature we find that the above equation gives 
rise to the continuity equation 

^ J eefj,\toi + Qoi)d^x = - j) [ee^*(tpi + ©pi)] dC', 

where the integration is on a three dimensional volume 
V bounded by the surface S. This leads to interpret tA 
as the energy-momentum tensor of the gravitational field 
[30, 32]. 

III. EXACT TELEPARALLEL TWO BODY 
PROBLEM 

The metric characterizing the space, with axial- 
symmetric static gravitational field in cylindrical coor¬ 
dinates (t, z, p, (fi) about the z axis can be written as 

{dz^ + dp^) — e~''p^dip^, ( 10 ) 
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where fi = fi{z,p), v = i'{z,p) and the coordinate while (13) gives rise to 
ranges are 

01^ = 0 => 

OO > t > —oo, OO > 2: > —00, 00 > p > 0, ZTT > (fi > 0. 




( 20 ) 


The metric (10) gives rise to the following diagonal tetrad 


e^=diag(e"/2, pe-"/^). (n) 


The field equations (7) read non-vanishing components 
of the energy-momentum tensor as 


2 2 r\ / \ p 

“t“ 2iyfXzz ^zz H" f^pp ^pp) 

^( 12 ) 

01^ = -02^ = (2pp -f 2z/p -h pz/2 - p;/2) ^ (13) 


00 ° = 




327r 


01^ = 02^ = 




IfiTTp 


{Pz +l^z- pt^zVp) : 


(14) 


03° 




(y^ + ‘i.h'zz + 2p zz H“ “1“ ^f^pp “1“ , 

(15) 


where pz = If: Mp = = C#; and so on. We 

solve the system for vacuum 0^-1 = 0. We combine 


00° - 03° = ^ {l^zz + Vpp + Vp/p) = 0. (16) 

OTT 

The vanishing of the quantity between the brackets in 
(16) is just the Laplace equation V'^v = 0 in cylindrical 
coordinates with an axis of symmetry about 2 . The field 
equation (16) gives the linear Newtonian gravity; then 
the field of two monopoles can be introduced as a su¬ 
perposition of the gravitational field of two singularities 
as 




2mi 2m2 
r_ r+ ’ 


(17) 


Combining (19) and (20), and substituting from (17) we 
evaluate 


A(2,p) 


mi 

ml' 

2 , mim2 

\(z^ + P^-d‘^\ 

- 1 

A 

A 



r_ 

+ 

1 



where the constants of integration have identified by im¬ 
posing the boundary condition; the space is asymptoti¬ 
cally fiat. Consequently, we have 


p(2,p) = 2 


mi m2 


'ml ml' 

- 1 - 

— 

A ' 4 

+ 

1 


r_ 


17111712 


r- r+ ) 


( 21 ) 


As mentioned above the contribution of the classical 
theory of gravity is represented by the linear equation 
(17), while the non-linear contribution is given by the 
quadratic terms in equation (21). The obtained solution 
is equivalent to the general relativistic one [2]. It has 
been pointed out that the space (10) violates the ele¬ 
mentary flatness in the region between the two masses 
[11, 33]. This can be seen when taking a circle with 
a center on the 2 -axis in the subspace 2 = constant, 
t = constant with | 2 | < d, one finds that the ratio of the 
circumference C to the radius R of that circle as i? —0 
is C/R 27re'^(^’°\ This requires that A(2,0) to van¬ 
ish everywhere on the plane which is not fulfilled for the 
2 -axis region with | 2 | < d, accordingly the ratio C/R 
does not approach 27r. This type of singularity is called 
conical singularity, the only physical interpretation from 
the GR point of view is to assume some sort of matter 
with Tpn ^ 0 within this region, i.e. to introduce a strut 
with zero active gravitational mass carrying the stress 
between the two masses [10, 11]. However, the TEGR 
solution might enable to examine the role of the telepar¬ 
allel gravity to reveal the nature of the repulsion between 
the two singularities. 


where nii and m 2 are two constants of integration, r_ 
and r_|_ are bipolar coordinates of an observer in {z,p) 
plane referred to two singularities, respectively, separated 
by a distance of 2d with the origin of a reference frame lies 
at the midpoint between the two singularities on the 2 - 
axis. The direction of the two singularities can be defined 
using the bipolar coordinates as 

= (2 ± d)^ + p^. (18) 

The linearization of the GR tells that goo is the Newto¬ 
nian potential; then the constants mi and m 2 represent 
masses. We introduce the function X = v + p, such that 
the field equation (14) reads 

01^ = 0 A 2 = pUpVz, (19) 


IV. VISUALIZATION OF BBHS SPACE-TIME 

In this section we investigate the singularities and hori¬ 
zons of the above solution via curvature verses torsion 
invariants. We evaluate the Ricci scalar of the BBHs 
space-time described by the induced metric (10) 

R{z,p) = -^^{‘^PPpp + 2ppzz + pn"/ -2np + pn'^). ( 22 ) 

It is clear that the solution given by (17) and (21) im¬ 
plies the vanishing of the Ricci scalar as expected for a 
vacuum solution. Recalling that the energy-momentum 
component (15) of vacuum and by comparison with (22), 
we find that the Ricci scalar is enforced to produce the 
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Laplace equation = 0 in vacuum. We conclude that 
the vanishing of Ricci scalar is natural for the vacuum 
choice. Consequently, all the invariants relevant to the 
Ricci would vanish and do not afford visual images of the 
space-time. Moreover, the solution is static, we expect 
all the magnetic part of Weyl tensor to vanish. This left 
us with only two possible non-vanishing invariants^, that 
are 


3-2/i 


WIR = - QpVpp{^,pVpp -I- - p^Vzp - Vp- Pp) 

+ - ^P^ZZ {(1 - PPp)Vp + PzVzP + Pp} + pIpP^ 

+ ‘^PPppi-vlp + PzzP - vIp + 2j^p) -f I2p^v% + pI^p^ 


pz 

- \2pvp^{-V:, - Pz + pi'zPp + PzVpp) -f y^p^ - 4,ypP 

- 2pp^z{ylp - ‘^Vp + vz^p) + v'^p^ + 3pl + 3/ip -I- dpz^z 
+ ^li'^^zP^ + 7 - QppP + iplp^ + iplp^) + yp {-6 pIp 

+ 6/ip — Ay'^p — 6plp) +y‘^{3 — QppP + 3/ipP^ -I- 3plp^)] , 


(23) 


and 

_ 

W2R = [pIpP^ + 2p(-r'pP -f pzzP - ylp + ‘^yp)Ppp 

- 9i^ppP^ + I8p(/ipi/pp -I- yzzP - PzyzP -yp- Pp)ypp 
d” PzzP ‘^Pi^pP 2iip -f y^p)pzz ^^zzP 
+ (18 {{-ppP + l)yp + PzyzP + Pp})pyzz - ^Gp^y% 
+ 36p{-y:z - Pz + pyzPp + Pzypp)ypz + y^p^ - ^ylp 
+ (-5 - %pIp^ + 2vIp^ - 9pIp‘^ + I8ppp)yl - 9pl 
+ {^8plp - 18/ip - -h 18plp)yp + y^p^ - 9pl 
+ (-9 -f 18/ipp - 9pIp^ - 9pIp^)vI - I8pzyz\ x 
{PppP +PzzP + ‘2.yp-vlp-ylp). (24) 

Substituting from (17) and (21) into (23) and (24), a 
straightforward calculations indicate that the the invari¬ 
ants Wli? and lT2i? diverge when r± 0. Since the 
solution has a curvature singularity at r± = 0 but goo 
does not vanish for r± > 0, the solution has no horizons. 
So the singularities at r± = 0 are physical and naked. 

On the other hand, the three independent Weitzenbock 
invariants h = T'^^^Tijk, h = T'-^^Tjik and h = 
where T* = might help to analyze the singularities. 
We omit using these invariants separately in this study 
due to the lack of invariance under LLT. However, their 
combination gives a well defined invariant, that is the 
teleparallel torsion scalar [34]. The contracted Bianchi 
identity (5) shows that if we want the teleparallel torsion 
scalar to vanish this implies the vanishing of the total 
derivative term as well. As a matter of fact, this cannot 
be happen, at least, for this solution. We show this in 


^ It is worth mentioning that the so called Kretschmann scalar is 
just 8WIR in this space [17]. 


(b) Torsion contour 

FIG. 1. Visualization of BBHs space-time: (a) a Holographic 
visualization of the teleparallel torsion scalar field formed by 
BBHs gravitational field; (b) The contour lines of the telepar¬ 
allel torsion show a repulsive zone in the inner region of the 
BBHs. The constants have been chosen as mi = 1, m 2 = 1 
and d = 2. 

more details, substituting from (11) into (1) and (4), we 
evaluate the teleparallel torsion scalar 

T{z, P) = ^ (2Mp + 2i^p - pvz - P^D ■ (25) 

Recalling that the energy-momentum component (13) 
of vacuum and by comparison with (25), the teleparallel 
torsion scalar of binary black holes reduces to 

T{z,p) = -yle-^^0. (26) 

The teleparallel torsion scalar vanishes if zzj, = 0. Conse¬ 
quently, the potential function v becomes a function of 
the coordinate p only, which contradicts the vacuum so¬ 
lution. So we conclude that the vacuum solution in fact 
prevents the torsion scalar to vanish completely unlike to 
the Ricci scalar. 

Substituting from (17) and (21) into (26), we find that 
the torsion scalar diverges when r± —> 0 which is in agree¬ 
ment with the Weyl invariants WIR and W2i? results. 
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However, the study via the teleparallel torsion invariant 
T is much simpler than W\R and W2R invariants. The 
non-vanishing torsion scalar of vacuum solutions, in gen¬ 
eral, might help to visualize the space-time properties. 
The case here of the BBHs vacuum solution can be seen 
in Figure 1. Recognizably, the contours of the torsion 
scalar similar to the magnetic field lines of binary pul¬ 
sars. However, the inner region shows clearly a repulsive 
behaviour of the torsion held lines. This recommends the 
non-vanishing teleparallel torsion scalar Held to interpret 
the compression force that keeps the two singularities 
apart. Finally, we Hnd that the singularity analysis of 
the BBHs via curvature verses torsion are in agreement 
as expected in the non charged or TEGR cases [34]. In 
the following section we use the acceleration tensor to 
investigate the repulsive gravity zones of the BBHs Held. 


V. REPULSIVE GRAVITY ZONES 


Consider the vierbein fields as reference frames of an 
observer in space-time defined by an arbitrary timelike 
worldling C. These frames can be characterized in an 
invariant way by antisymmetric acceleration tensor 
This tensor generalizes the inertial accelerations of the 
frame, in analogy with Farady tensor, the acceleration 
tensor can identify —>• (a, Jl), where a is the transla¬ 
tional acceleration = “(m)) ^ frequency 

of rotation of the local spatial frame w.r.t. non-rotating 
frame [35]. 

Let the worldline C of the observer be denoted by 
x^{t), where r is the proper time of the observer, whose 
frame is adapted such that its components identify the 
velocity and acceleration along C respectively by e(o) * = 
M* = dx^/dr and a® = Du^/dr. The absolute derivative 
D/dr is w.r.t. the Christoffel symbols F®;. as 


a 


I 


Dvd 

dr 


= U'^Vje(o) 


d?x^ " j dx^ dx’^ 
dr‘^ dr dr 


(27) 


We see that if u® = e(o) ® represents a geodesic trajectory, 
then the frame is in free fall and a® = 0 = 4>(s>){^i)- There¬ 
fore we conclude that non-vanishing values of the latter 
quantities represent inertial accelerations of the frame. 


a = 


Du^ 

dr 


De 


(0) 


dr 


= M-^Vo-e 


je(o) 


(28) 


Assuming that the observer carries an orthonormal 
tetrad frame ®, the acceleration of the frame along the 
path is given by 


= (29) 

where is the antisymmetric acceleration tensor. It 
follows from Eq. (29) that 


Dcu 


= e 


dr 


= e" i M-’V,e^ 


(30) 



(a) ^(0)(0)(l) component 



(b) T(o)(o)( 2) component 


FIG. 2. acceleration tensor components: (a) along the 2 - 
direction have some regions with negative values near the sin¬ 
gularities and between the two singularities; (b) along the p- 
direction have only positive values. The constants have been 
chosen as mi = 1 , m 2 = 1 and d = 2. 


The acceleration vector a® defined by Eq. (28) may be 
projected on a frame yielding 

= u i a® = i Vje(o) * = ^!>(o) ^ ■ 

It is possible to rewrite the acceleration tensor in the 
form 

+ ^iJ.{0)v — Ti,(o)ij] 1 ( 31 ) 

where T^up — and djCpi is 

the torsion tensor of the Weitzenbock space-time. The 
expression of (jjpi, is invariant under coordinate transfor¬ 
mations [35]. The values of the 6 components of the ac¬ 
celeration tensor may be used to characterize the frame, 
since (ppi, is not invariant under local SO(3,l) (Lorentz) 
transformations. Alternatively, the frame may be char¬ 
acterized (i) by the identification u® = e(o) ® (this equa¬ 
tion fixes 3 components, because e(o) ^ is fixed by nor¬ 
malization), and (ii) by the 3 orientations in the three- 
dimensional space of the components e(i)®, 6 ( 2 )®, 6 ( 3 )®. 

Eor a frame that undergoes the usual translational 
and/or rotational accelerations in flat space-time, (jjpi, 
yields the expected, ordinary values. An interesting ap¬ 
plication of the acceleration tensor is the following. 
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Let us consider a static observer in the BBHs space- 
time described in the cylindrical coordinates (t, z, p, ip) by 
the vierbein (11). The frame of the observer must sat¬ 
isfy 6 ( 0 ) * = 0 = M*, the translation accelerations identify 
acceleration tensor components 

<('(0)(/.) = T’(0)(0)(/x) = e(0)*e(^)^r(0)ij, 

where the non-vanishing components of the translation 
acceleration tensor of the BBHs space-time are 

The translation inertial acceleration plots, Figure 2, in¬ 
dicate that there are repulsive gravity regions along z- 
axis which is identified by 0(o)(i) < 0 values, while they 
indicate that the gravitational acceleration is always at¬ 
tractive along p-axis. In this case the translational ac¬ 
celeration vector a = (<(>(o)(i), </i’(o)( 2 )i 0; O) can be given 

by 



FIG. 3. The acceleration in the BBHs gravitational field 
shows a deceleration region between the two singularities. 
This indicates a repulsive gravity region. The constants have 
been chosen as mi = 1, m 2 = 1 and d — 2. 



e p/'^z + -iype 


-A‘/ 2 ; 


where z and p are unit vectors along z and p directions, 
respectively. Consequently, 

/ 2 2 \ 1 /^ 

+ (32) 

For simplicity we take m-i ~ m 2 , along the z-axis outside 
the BBHs, i.e. |z| ^ d. The acceleration is 


^^out — 2 


(mi-bm2)zd ^^(,,.o)/2 
(z + d)2(z-d)2 


g Wim 2 


Also, its value in the inter-medium region, i.e. |z| d, 
is given by 


^ (mi-bm2)zd „_y(^,o)/2 
(z + d) 2 (z-d) 2 ® 

^ ^ rnim2 


The above expressions represent the inertial acceleration 
needed to cancel exactly the gravitational acceleration 
maintaining the frame static. It is clear that the outer 
regions of the BBHs along z-axis experience the usual 
attractive gravity, while the inner region indicates a re¬ 
pulsive gravity zone. We examine the acceleration along 
p-axis, we get 


also interested to examine the acceleration at the center 
of mass, which can be given as 

a(0,0)= 

The inertial acceleration at the center vanishes if mi = 
m 2 . This is expected as the gravitational acceleration 
vanishes at the same condition so that the inertial ac¬ 
celeration must vanish to maintain the frame static. In 
Figure 3, we plot the acceleration in (z,p) plane of the 
BBHs. The plot shows that acceleration increases as an 
observer falls towards any of the two monopoles, while it 
decreases at some regions between them. This deceler¬ 
ation regions identify the repulsive gravity zones of the 
BBHs space-time. 


VI. CONSERVED QUANTITIES OF THE 
BINARY SOLUTION 

The energy localization problem has been studies 
within the tetrad formalism in many literature [36, 37]. 
The Lagrangian of the TEGR using forms techniques 
takes the form^ 

V = -Y^r“ A* - 2 1 , ( 33 ) 

where + denotes the Hodge duality in the metric gap 
which is assumed to be flat Minkowski metric pap = 


a(0,p) 


(mi +m 2 )p p(o^p )/2 

(p2+d2)3/2 


mi -I- m 2 


P 


2 


The last term is evaluated where p^ d, which indicates 
that the gravity along the p-axis is always attractive. We 


® The role of the addition of non- Riemannian parity odd pseu- 
doscalar curvature to the Hilbert-Einstein-Cartan scalar curva¬ 
ture has been investigated (cf., [38—40] and references therein). 
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Oap = rf*a5(+l, —1, —1, —1), that is used to raise and 
lower local frame (Greek) indices. The variation of the 
total action w.r.t. the coframe gives the equation of mo¬ 
tions in the from [41] 

- fa = Sc, where 

0'0°‘ 


that is the canonical energy-momentum current 3-form 
of matter which is the source. Respect to the general 
Lagrange-Noether scheme [42], one can derive from (33) 
the translational momentum 2-form and the canonical 
energy-momentum 3-form: 


■“ dT^ ~ 8 ^ * 
dV 

■“ ~ CaJ V -f (eaJT^) A Tip. 


(35) 

(36) 


where T^ is the purely Riemannian connection and 
is the contorsion 1-form that is connected to the torsion 
by the relation 


r“ ‘^=- /C“3 A 1?/^. 


(41) 


The teleparallel model (33) belongs to the class of quasi¬ 
invariant theories. One can easily show that under a 

change of the coframe = ^A(x)“^^ the La- 
grangian (33) changes by a total derivative: 


V(d') = V(d) - —d 




(42) 


Besides to (42), it is an easy task to check that Eq. (40) 
varies like 


Because of the geometric identities [43], the Lagrangian 
(33) can be rewritten in the form 

V = -ir“A7^„. (37) 

The existence of the connection field T°‘p plays an im¬ 
portant role in the regularizing due to the following: 




- ^—d 

IOtt 

£\{d') = (A-i) 
~ 1 ^^ 


^j0{d)+d{A-^f^An0{^) 

\A-^)\{A-^)\dA\Ar^p\ 


.(43) 


i) The TEGR theory be invariant under the LET of the 
coframe, i.e., the Lagrangian (33) is covariant under the 
change of variables 

d'“ = A°‘pd^, 

T'J = A^^J:/{A-y^ - {A-^)\dA^^. (38) 

Because of the non-invariant transformation law of Tq,^ 
as Eq. (38) shows it will not vanish in any other frame 
connected to the first by a LLT. 

ii) The connection plays an essential task in the 
teleparallel framework. This task describes the inertial 
causes which occur from the selection of the reference 
system [41]. The evolvements of this inertial in many 
situations lead to non-physical results of the total energy 
of the system. Therefore, the function of the teleparallel 
connection, is to deviate the inertial involvement from 
the really gravitational one. Because of the teleparal¬ 
lel curvature is vanishing identically, the connection is a 
“pure gauge”, that is 

= {A-y^dA^^. (39) 

The Weitzenbdck connection forever has the form (39). 
The translational momentum of Lagrangian (33) has the 
form [41] 

'Ha = ^ lap-y, 

Ta^ djf- fd _ ,Ca^ (40) 


The total conserved charge in the TEGR theory has 
the form [43] 


Jit^) 


1 

IGtt 



1 

IOtt 


f (44) 

JdS 


with Using LLT Eq. (44) changes as 


J'{.L A lla\ 

with = (^A~^Y^dA ^The vector ^ = C,^di is inde¬ 
pendent of the choice of the frame while its components, 
i.e., change as a vector. Now let use apply Eq. (44) 
to the tetrad (11). Using the spherical local coordinates 
(t, r, 0, (fi) the binary solution, by proceeding the follow¬ 
ing transformation {z = rcosd, p = rsind). Recalling 
the vierbein (11), the frame is described by the coframe 
components: 


= e 2 dt, 

(46) 

= e'‘‘2' ’dr, 

(47) 

= re 2 d9, 

(48) 

Q —u(r,6) 

di =rsinde 2 dip. 

(49) 


If we take coframe (46) we get the non-vanishing compo- 
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-1 V,. u{T,a)-p.{r,e) . 

Tn = -5- dt. 


-L n — 


' dt. 


nents of the Riemannian connection f £ in the form 

(50) 

(51) 

(52) 

(53) 

'dLp. (54) 


~2 dedr — 2rd9 - r^fij-dO 

i 1 = - 


2r 


~3 , {2-rVr) 

r^ = —--- -e 2 r dtp, 


3 , ^2c0t9—l'g -./(r,8)-^(r,fl) 


r 2 = — sin0 


Using (50) as well as the Riemannian connection and 
substitute into (40) we finally get 




Hq = sin9e ^ < r[4: + r{fj,r — Vr)]d9 A d(p 


— [2 cot 9 — vg + pgjdr A d(p 


~ -M(r,e) 

'Hj=2cos0e 2 dtAdif, 

~pL(r,6) 

'H2 = 2sin0e = dipAdt, 

u(r,a) r[2 + rvr + rp,r]d9 Adt + [vg + d,g]dt A dr 

iLo — 6 ^ • 

J r 

(55) 

For the BBHs solution we have for a vector field ^ = 
with constant holonomic components, C®, 
in the coordinate system used in (46), the direct evalua¬ 
tion of the integral (44) yields the following infinite total 
conserved charge 

J = f CHa = {mi +m2-r) C° 

lOTT Jqs 


3r 

4 f m^ — m\ + 2mim2 ^ 


-) 


C“ + ohr (“) 


One possible solution of the above unfamiliar result is the 
use of regularized form of total charges 

1 


j(e,^) = 


1071 


dS 


^p/37 


^fmi=0, m2=0]) ^ da/37- 


(57) 

Using (57) we get the total conserved quantities in the 
form 

J = = {mi + m 2 ) C° 

lOTT JdS 

' m| — m^ — 2 toiTO 2 \ *0 

3 ^^ ) ^ 

which is a consistent with the results in [44]. It is 
clear that the leading term is consistent with the case 
of Schwarzschild field, which presents the total mass of 
the system mi + m 2 - 


The analysis of the total conserved quantities are com¬ 
bined only with the C*^-component, while other terms 
vanish. Consequently, we conclude that the total con¬ 
served quantities are purely an energy. Since we have as¬ 
sumed that the separation between the two point masses 
is constant, we do not expect any energy loss in the orbit. 
But in order to estimate an initial guess for the energy 
loss by a binary system due to its gravitational radiation. 
We consider the case mi = m 2 = m and the separation 
between the two point masses to be a function of time, 
i.e. d = d{t). Using (58) the energy contained by the 
system E{t) can be written as 


E{t) 


2 m^ 

2m- - - 

3 r 


8 m‘^d{t)'^ 

15 


(59) 


The rate of the energy loss by the binary system is 


E 


16 m? d 
15 


d. 


Using (59) and (60) we evaluate 


(60) 


E 

E 



(61) 


which agrees with the Keplerian treatment of orbital pe¬ 
riod decay [45]. It is worth to mention that the solution 
is an exact one for which we can find such quantities of 
physical interest as radiation patterns without assuming 
a weak field initial data. 


VII. CONCLUSION 

In the frame of the TEGR theory, we have derived an 
exact solution of BBHs space-time. The field equations 
enforce the Ricci scalar to vanish, while they prevent the 
vanishing of the torsion scalar. This property has been 
used to visualize the space-time of the BBHs field. This 
fact can serve in the torsion vs curvature singularity anal¬ 
ysis. As a matter of fact, in the Riemannian picture, all 
the invariants due to the Ricci tensor vanish, while the 
only non-vanishing components of the Weyl (dual) tensor 
are the electric components. So we used the invariants 
WIR and W2R in this study. Additionally, we used the 
teleparallel scalar, in the Weitzenbbck picture, to ana¬ 
lyze the singularities. The results identify the physical 
singularities at r± = 0. However, in our case here, the 
torsion contours appear as repulsive lines, which might 
afford a physical explanation to the stress keeping the 
system balanced. 

We have calculated the acceleration tensor to study 
the inertial acceleration required to keep the frame in 
a certain inertial state in the BBHs space-time so that 
in static frames the the inertial acceleration is precisely 
opposing the gravitational acceleration. We have identi¬ 
fied repulsive gravity zones via the acceleration tensor in 
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the inter-medium keeping the system balanced. In this 
way, we do not need to introduce a strut or membrane 
to interpret the configuration as in the GR theory. 

We have studied the conserved quantities of the BBHs 
gravitational field. The regularized total charge calcula¬ 
tions have been used to estimate the rate of the energy 
loss due to emission of gravitational radiation by the sys¬ 
tem. 

In conclusion, we would like to point out that the non¬ 
vanishing torsion scalar of the vacuum solution of BBHs 
space-time can be used to study the singularities of this 
solution. Also, the modification of the TEGR Lagrangian 
by adding a scalar field Lagrangian. In some cases the 
scalar field can be induced directly from the torsion [46] , 
which can be used to interpret the repulsive gravity zones 
near the naked singularities and in the inter-medium of 
the BBHs space-time. Actually, in this paper we just 
alarm the importance of studying these types of solu¬ 
tions in the TEGR framework. However, the singularity 
and the gravitational energy-momentum tensor analysis 
in the TEGR framework still needs more investigations. 
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Appendix A: Notation 


As is known, the exterior product is indicated by A, 
however the interior product of a vector ^ and a p-form T 
is indicated by ^J'l'. The vector basis dual to the frame 
1-forms ■(?“ is indicated by Cq, and they fulfil = 

5^. By using local coordinates a;*, we get = hfdx^ 
and Ca = h'^adi with hf and are the covariant and 
contravariant components of the tetrad field. The volume 
4-form is defined by rj A A A . Additionally, 

by using the interior product we define 


def. I 

Va = ea\r] = 


1 


ea/3-/5A A 


with Cap-fS is completely antisymmetric tensor and 

£0123 = 1 - 

l/a/3 ■— — '^£a/37(5 A d , (-^1) 

l?a/37 ■— £ 7 ] d , (-^^) 

that are the bases for 3-, 2- and 1-forms respectively. 
Einally, 


XjaPfiu — — Ci/J £/3j 

is the Levi-Givita tensor density. The ry-forms fulfil the 
following useful identities: 

d^ A 77a := df 77 , 
d^ A 77;,^ := S^r]^ - 
d^ A 77a/ii/ -f S^TJua T d^77a/i, 

d^ A := . (A3) 

The line element ds^ := gap'd°‘ ^ d^ is fulfil by the space- 
time metric g^p- 

One can consider teleparallel geometry as a gauge the¬ 
ory of translation [42, 43, 47-49]. In such geometry the 
coframe d“ plays the role of the gauge translational po¬ 
tential of the gravitational field. General relativity can 
be reconstructed as the teleparallel theory. Erom geo¬ 
metric viewpoint, teleparallel gravity can be regarded 
as a particular case of the metric-affine gravity whose 
coframe 1-form d“ and Lorentz connection are due to 
distant parallelism constraint Ra^ = 0 [43, 49, 50]. In 
such geometry the torsion 2-form 

r“ = i:)d“ = dd“ -b Tp’^ A dd 

= A d" = A da:^ (A4) 

occurs as the gravitational gauge field strength, Ta^ is 
the Weitzenbock connection 1-form, d being the exterior 
derivative and finally D is the exterior covariant deriva¬ 
tive. The torsion can be divided into three irreducible 
parts: the tensor part, the trace, and the axial trace, pro¬ 
vided by [41, 47, 48, 51-53] 

W'J-a -y-a _ (2)-y-a _ {3)qpa ^ 


with 


(2)r“ := ^7?“ A r, (A6) 

o 

where T = {ep\T^) and CaJT = vectors of torsion 
trace, and 


(3)r“ := 


(A7) 


with V = A Tp) and = T^'^^'q^vXa the axial 

torsion trace. 
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